Time-resolved spectra of polar-polarizable chromophores in solution 
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A recently proposed model for steady-state spectra of polar-polarizable chromophores is extended 
to describe time-resolved spectra. The model, based on a two-state picture for the solute and 
on a continuum overdamped description for the (polar) solvent, grasps the essential physics of 
solvation dynamics, as demonstrated by the comparison with experimental spectra. The solute 
(hyper)polarizability is responsible for spectroscopic features that cannot be rationalized within 
the standard picture based on a linear perturbative treatment of the solute-solvent interaction. 
In particular, the temporal evolution of band-shapes and the appearance of temporary isosbestic 
points, two common puzzling features of observed spectra, are natural consequences of the molecular 
hyperpolarizability and of the consequent coupling between solvation and vibrational degrees of 
freedom. 



I. INTRODUCTION 



Time-resolved spectroscopy is a powerful tool to inves- 
tigate solvent dynamicsii^, one of the central issues in 
modern chemical physics, and a key to understand the 
rate of chemical reactions in solution^'"*. Detailed infor- 
mation on solvent dynamics can in fact be obtained from 
time-resolved spectra of properly chosen probe molecules 
in solution. Good solvation probes have strongly sol- 
vatochromic and intense absorption (and/or emission) 
bands in the visible region^, as not to interfere with tran- 
sitions in the ultraviolet region. Polyconjugated chro- 
mophores with good electron-donor (D) and acceptor (A) 
groups (also called push-pull chromophores) work fine in 
this respect: intense and low-lying transitions are guar- 
anteed by polyconjugation, and the large charge displace- 
ment from D to A (or viceversa) driven by the transition 
is responsible for a large variation of the dipole moment 
upon excitation, that makes the frequency of the transi- 
tion largely dependent on the polarity of the surrounding 
medium. 

In the simplest model for solvation, the solvent is de- 
scribed as a dielectric continuum medium hosting so- 
lute molecules into Onsager-type cavities^. The electric 
dipole moment of the (polar) solute polarizes the sol- 
vent so that the solute itself experiences an electric field, 
the reaction field, proportional to the solute dipole mo- 
ment via a proportionality constant that depends on the 
macroscopic properties of the solvent, and on the shape 
and size of the Onsager cavity. Two contributions to 
the reaction field can be distinguishedSi^. The electronic 
clouds of solvent molecules surrounding the solute de- 
formate in response to the solute dipole moment. The 
resulting electronic polarization is very fast with respect 
to the relevant degrees of freedom of the solute. As a con- 
sequence, the dynamics of the electronic polarization is 
irrelevant to our discussion^. The second contribution to 
the reaction field originates from the reorientation of the 
(polar) solvent molecules around the solute. This contri- 
bution, that of course vanishes in non-polar solvents, is 
slower than the solute degrees of freedomSii. It is respon- 
sible for the large solvatochromic behavior of push-pull 



chromophores^ and dominates the temporal evolution of 
observed spectra in time-resolved experimentsi. 

In a typical time-resolved experiment the system is 
shot by a pulse, or a series of pulses, that induces a transi- 
tion between the ground and the excited state (or vicev- 
ersa). During this process the solvent molecules have 
not time to reorient in response to the sudden change 
of the solute dipole moment, and an out of equilibrium 
(orientational) polarization is generated. Provided the 
state generated by impulsive excitation lives long enough, 
the orientational contribution to the reaction field {For) 
slowly evolves from its initial value (i.e. that equilibrated 
with the solute dipole moment before excitation), to- 
wards the equilibrium value for the state reached upon 
excitation. Within linear perturbation theory, the energy 
of the states involved in the transition is linearly affected 
by the reaction field, and the frequency of the relevant 
transition acquires a linear dependence on the reaction 
field: hijj{For) = hLu{For = 0) — For{l^E — Mg)j where 
fJ-GiE) is the electric dipole moment of the ground (ex- 
cited) stated. Then, the experimentally accessible tem- 
poral evolution of the transition frequency, u}{t), gives 
direct information on the temporal evolution of the ori- 
entational solvent polarization. 

The continuum dielectric model reduces the highly ir- 
regular and complex motion of the solvent molecules to 
the motion of a single effective solvation coordinate, the 
reaction fieldiS. The dynamical behavior of this special 
coordinate is subject to large frictional random forces, 
that guarantee for an energy dissipation mechanism^. In 
the hypothesis of large friction (overdamped motion) the 
resulting dynamics is very simple. The solvation corre- 
lation function (also called the spectral correlation func- 
tion). 



C{t) 



Uj(t) — Uj{(X)) 

uj{0) — uj{co) 



(1) 



has in this hypothesis a simple exponential form, and the 
corresponding relaxation time, ts, coincides with the lon- 
gitudinal relaxation time, r^, of the pure solvent^. Devi- 
ations from this simple behavior are expected in the case 
of low- friction or when the translational motion of solvent 
molecules dominates over orientational motioniii. 
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Extensive experimental work has been devoted by sev- 
eral groups to collect time-resolved spectra of polar chro- 
mophores in solution^^i^'^ii^ii.'^iifiiiTii^i ^",2^22,28^ The 

interpretation of experimental data proved more diffi- 
cult than described in the simple picture summarized 
above. First of all a finite contribution to the Stokes 
shift (a;(0) — w(oo)) is given by the vibrational degrees 
of freedom of the solute molecule, as demonstrated most 
clearly by the observation of finite Stokes shifts for chro- 
mophores dissolved in non-polar solvents. Internal vi- 
brational coordinates are in fact slower than electronic 
degrees of freedom, and hence contribute to the Stokes 
shift. However, thanks to the different timescales of vi- 
brational and solvation motions, their contributions to 
the temporal evolution of spectra can be at least approx- 
imately separated^^. The basic idea is that, after im- 
pulsive excitation, vibrational degrees of freedom relax 
almost instantaneously with respect to solvation time- 
scales, so that vibrational relaxation is essentially done 
before solvent relaxation starts. Accordingly, the zero- 
time frequency, w(0), to be inserted into Eq. is not 
the true zero-time frequency, i.e. that relevant to the un- 
relaxed (vertical) system, but represents the frequency 
measured at an effective zero-time for the solvation mo- 
tion, i.e. the hypothetical time when vibrational coordi- 
nates are fully relaxed, but the solvent is still frozen in 
its vertical configuration. This effective t — time, first 
introduced by Maroncelli and coworkers in the analysis of 
time-resolved fluorescence spectraii, has no experimen- 
tal counterpart, and a;(0) cannot be measured directly. 
Maroncelli proposed a clever approach to extract w(0) 
from experimental steady-state spectra collected in polar 
and apolar solvents^*. The procedure unavoidably intro- 
duces uncertainties in the estimated w(0) that add to the 
uncertainties in the frequencies estimated from the broad 
and asymmetric fluorescence bandsi^. Precise estimates 
of C{t) are therefore difficult. On a more fundamental 
vein, the simplest solvation model outlined above, and, a 
fortiori, the approach suggested by Maroncelli, both rely 
on the hypothesis of instantaneous vibrational relaxation 
and on the hypothesis of an ideal solvation probe: i.e. of 
a molecule whose spectral properties are largely affected 
by the surrounding medium, but whose presence does not 
affect the solvent dynamics. The assumption of a com- 
plete separation of vibrational and solvation time-scales 
is expected to work well unless one is interested in the 
very early times of solvation (the first few tens of fem- 
toseconds). Deviations can be predicted for molecules 
with very slow (conformational) degrees of freedom, par- 
ticularly if solvents with very fast relaxation dynamics 
are considered. The assumption of the solvent response 
being unaffected by the perturbing solute is only valid 
in the linear regime for the solute-solvent perturbation. 
But linear approaches are inadequate to treat solvation 
of push-pull chromophores. Just in view of the strong 
absorption in the visible region, in fact, push-pull chro- 
mophores have large polarizabilities (and hyperpolariz- 
abilities): the reaction field then not only affects the 



solute energy, but also its dipole moment. The dipole 
moment in turn affects the reaction field in a feed-back 
mechanism that is responsible for large non-linearity. 

In recent years we have developed a simple yet power- 
ful method to describe steady-state electronic and vibra- 
tional spectra of polar chromophores in solution2SiSSiSL2£. 
The solute is described in terms of two electronic states 
linearly coupled to internal vibrations (Holstein coupling) 
and to the reaction fieldS. This model allows for non- 
perturbative solutions and naturally describes several 
features of experimental spectra that were not under- 
stood in standard linear approaches. The reliability of 
the model is proved by extensive comparison with steady- 
state electronic and vibrational spectra collected for sev- 
eral chromophores in different solvents^SiSSi^. In this pa- 
per we extend the model to describe time-resolved spec- 
tra. In particular we show that accurate time-resolved 
spectra of polar-polarizable chromophores can be cal- 
culated adopting the same set of parameters used for 
steady-state spectra without the need to introduce ad- 
ditional adjustable parameters. Working in the hypoth- 
esis of a complete separation of solvation and vibra- 
tional dynamics, the model accounts for deviations of 
C{t) from the simple exponential form, and rationalizes 
the observation of time-dependent band-shapes and the 
appearance of temporary isosbestic points (TIPs) . These 
characteristic features of time-resolved spectra of polar- 
polarizable chromophoresi2i22i2ii22iS^, are a direct conse- 
quence of the solute (hyper)polarizability, and cannot be 
rationalized within the standard linear picture for sol- 
vation outlined above. As a matter of fact these fea- 
tures have often been taken as evidences of the failure 
of the continuum solvation modeliii^ii^ or for the fail- 
ure of the two-state description of the chromophorc^'^. 
We demonstrate instead that they can appear as natural 
consequences of the molecular polarizability and can be 
described within the simplest model for solvation, pro- 
vided the fundamental interactions of electronic degrees 
of freedom with molecular vibrations and with the sol- 
vent reaction field are treated beyond linear perturbative 
approaches. 



II. THE MODEL 

We consider two electronic states linearly coupled to 
molecular vibrations and to solvent degrees of freedom, 
according to the following HamiltonianS: 



n 



2zop - t ^ 



'■Dj 



E 



,2,o2\ 



(2) 



where the first two terms describe the electronic prob- 
lem, with p — a\ ijClam counting electrons on A site, 
creating an electron with spin tr on D 



and a 



A, IT' "-0,(7 



and A site, respectively. Qi and Pi, with z > 0, are 
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the coordinates and conjugated momenta describing the 
coupled vibrations of the chromophore, and Qq is the 
effective solvation coordinate. All coordinates are har- 
monic, with frequency uji] the relevant relaxation energy 
is In the following we explicitly account for a sin- 
gle coupled vibration, Qi. This simplifies equations and 
limits the number of adjustable parameters in the model, 
with no major loss of accuracy in the calculation of elec- 
tronic spectra™. The extension to the multimode case 
is possible, and interesting effects of the mixing among 
different vibrational coordinates are indeed expected in 
vibrational spectrEkSL^. 

In the adiabatic approximation, the Hamiltonian in 
Eq. can be diagonalized at fixed Q^'s on the basis 
of \DA) and \D+A^) states^i^. The resulting ground 
and excited states are fully defined in terms of a single 
Qi-dependent parameter, p — {G\p\G), as follows: 



\G) 
\E) 



./l^p\DA)+^\D+A-) 
-^\DA) + ^l~p\D+A-) 



(3) 



The corresponding energies define the multidimensional 
ground and excited state potential energy surfaces (PES), 
as shown in Fig. la. The calculated PES are interest- 
ing in several respects. All coordinates inserted into the 
Hamiltonian are harmonic, and no term appears into 
Eq. (0) explicitly accounting for the mixing of different 
coordinates. The PES in Fig. 1 are instead clearly an- 
harmonic and indicate a large mixing between the two 
involved coordinates. The softening of the ground state 
PES, and the hardening of the excited state PES are also 
recognized. All these effects are due to the coupling of 
the Qi's to the electronic systemSi. In particular we un- 
derline that the coupling between vibrational and solva- 
tion coordinates stems out naturally from their common 
coupling to the electronic system. 

The PES in Fig. 1 represent the potential for the mo- 
tion of slow coordinates, Qi's. Internal vibrations with 
typical frequencies in the mid-infrared region, are true 
quantum-mechanical objects. The relevant anharmonic 
problem defined by the PES in Fig. 1 has been solved nu- 
merically, and large effects of the anharmonicity of inter- 
nal vibrational modes have been demonstrated in static 
hyperpolarizabilities^SiS. On the contrary, vertical pro- 
cesses, like, e.g., absorption and emission are hardly af- 
fected by anharmonicit}i22i24i2^. In particular, vertical 
processes are quite accurately described by replacing the 
two anharmonic PES by the two parabolas that best fit 
the ground and excited state PES in the relevant region 
of the Q space (best harmonic approximation22iSiM). Of 
course the role of anharmonicity can be recognized in the 
variation of the effective harmonic parameters fitting ab- 
sorption and emission processes in this local-harmonic 
approximation. 

Within the local harmonic approximation, electronic 
spectra and their variation in the Q-space are fully de- 
fined in terms of p. The vertical transition frequency and 



the transition dipole moment ar o^^i^^ : 



CT 



fJ-CT 



(4) 



where /ig is the dipole moment of ). The Huang- 

Rhys (HR) factor for the coupled vibrational coordinate 
(Qi) is2^: 



Ai^^(l-2pf 
Tiuji 



(5) 



Based on Eqs. I@J and ISJ, absorption and emission spec- 
tra at any specified position in the Q space are calculated 



as2i: 



— \/n'. 



(6) 



where w„ is the frequency of the n-th vibronic line, and 
I{uj) is the intrinsic line-shape function that we choose 
as a Gaussian with standard deviation a. M = 1,3 for 
absorption and emission, respectively. 

The solvation coordinate describes the slow orienta- 
tional motion of polar solvent molecules around the polar 
solute. The relevant coordinate, Qo, can then be treated 
as a classical coordinate, and, at finite temperature, it 
is responsible for inhomogeneous broadening. Since the 
solvation motion is coupled to electronic degrees of free- 
dom, and, through them, to internal vibrations, inhomo- 
geneous broadening has in these systems quite peculiar 
features. Due to thermal disorder in fact, the solution 
can be thought of as a collection of chromophores each 
one surrounded by a slightly different Qo configuration. 
Both electronic and vibrational motions are very fast if 
compared to solvation, so that each solute molecule read- 
justs instantaneously its polarity and its geometry in re- 
sponse to the local Qq configurations!. Thermal disorder 
in Qo is then responsible for disorder in electronic and 
vibrational properties, that shows up in electronic and 
vibrational spectra, and, most apparently, in the anoma- 
lous dispersion of resonant Raman frequencies with the 
excitation lineSLSS. With reference to Fig. 1, the ther- 
mally equilibrated state in either the ground or excited 
state PES is defined in terms of a Boltzmann distribu- 
tion of states lying along the dashed lines drawn across 
the isopotential lines in Fig. lb and Ic to mark the local 
minimum path with respect to the vibrational coordi- 
nate. The steady-state absorption and emission spectra 
are then sums of the contributions in Eq. © calculated 
for each point along the minimum path in the ground 
state and excited state PES, respectively, and weighted 
by the relevant Boltzmann distribution2L2&. 

This simple model has been applied successfully to de- 
scribe electronic and vibrational spectra of several chro- 
mophores, reproducing many features of experimental 
spectra that were not understood before22i22i2&. Just as 
an example. Fig. 2A and 2B shows the absorption and 
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emission spectra calculated in the proposed approach to 
fit experimental spectra of DCM dissolved in several sol- 
vents. These spectra accurately reproduce experimental 
data^S., and deserve some comment. The red-shift of ab- 
sorption and emission bands with increasing solvent po- 
larity (positive or normal solvatochromism) is typical of 
chromophores with a neutral ground state (p < 0.5)^. 
For these chromophores the solute polarity increases 
upon excitation, and the excited state is stabilized by 
a polar surrounding medium. The evolution of the band- 
shape with the solvent polarity is a more subtle phe- 
nomenon that can be understood only if the polarizabil- 
ity of the solute is accounted fori2L2&. In fact, the polarity 
of the (polarizable) solute increases with the polarity of 
the surrounding medium. For neutral chromophores then 
p increases towards 0.5, and HR factors decrease: even 
if inhomogeneous broadening increases with the solvent 
polarity, absorption and emission bands apparently nar- 
rows, since the underlying vibronic profile narrows (Ai 
decreases, see Eq. Q). For zwitterionic chromophores 
{p > 0.5), instead, the ionicity increases towards 1, and 
HR factors increase with the solvent polarity: in this case 
the effects of solvent polarity on inhomogeneous broad- 
ening and HR factors sum up to widen absorption and 
emission bands. This behavior is experimentally verified 
in dyes with inverse solvatochromismSUiSi^, and is re- 
produced in Fig. 2C and 2D for a hypothetical dye with 
the same molecular parameters as DCM, but a different 
zq, as to have a zwitterionic ground state. 

Another direct consequence of the molecular polariz- 
ability is recognized in the non-specularity of steady- 
state absorption and emission bands^^. The behav- 
ior is similar for neutral and zwitterionic chromophores, 
with narrower emission than absorption bands. In all 
chromophores, in fact, the mixing between the two basis 
states increases in the relaxed excited state, so that, in 
both cases steady-state emission spectra correspond to 
states with p nearer to 0.5 than states relevant to ab- 
sorption spectra. According to Eq. (jsj, HR factors for 
emission are always smaller than for absorption then jus- 
tifying quite naturally the observed narrowing of emission 
bands. 



III. RELAXATION DYNAMICS 

The success of the two-state model in reproducing 
steady-state absorption and emission spectra invites 
modeling time-resolved experiments, as related to the 
dynamics of non-equilibrium states created by impulsive 
excitation in either the ground or excited state PES. We 
are not interested in the fast relaxation of vibrational 
coordinates (typically accomplished within few tens of 
femtoseconds after excitation)^^*^^, rather we model the 
slower solvation dynamics in the hundreds of femtosec- 
onds to picosecond temporal window. The problem is 
therefore greatly simplified if we work in the hypothesis 
of a complete separation of vibrational and solvation dy- 



namics. First of all, just after a (sequence of) pulse(s) has 
created an out of equilibrium state in one of the PES, an 
instantaneous relaxation of vibrational coordinates takes 
place before solvation relaxation starts^^. Moreover, dur- 
ing the subsequent relaxation of solvation degrees of free- 
dom, internal coordinates immediately readjust with re- 
spect to the instantaneous configuration of the solvent. 
Due to the non-linearity of the coupled problem, the elec- 
tronic wavefunction readjusts following Qo and, as a con- 
sequence, the equilibrium geometry of the chromophore 
varies with Qq. Since electronic and vibrational degrees 
of freedom are much faster than solvation, they are al- 
ways in equilibrium with the instantaneous configuration 
of the solvent, and the relaxation of the effective solva- 
tion coordinate does not proceeds along a horizontal line 
(fixed Qi) in Fig. lb or Ic, but follows the dashed lines 
drawn across the isopotential lines to mark the minimum 
energy path in each surface. 

In a time- resolved fluorescence experiment, an impul- 
sive excitation drives the system from the thermally equi- 
librated ground state (centered on point G in Fig. lb) 
to the vertical excited state (point G' in Fig. Ic). Af- 
ter excitation, the instantaneous relaxation of vibrational 
degrees of freedom drives the system to point F, i.e. 
to the local minimum in the excited state PES, at the 
Qo relevant to the equilibrium ground state configura- 
tion. During the G— i'G'^F path, the Qq configuration 
is unaffected. Once the system has reached point F, sol- 
vent relaxation starts. Since vibrational degrees of free- 
dom follow adiabatically the solvent relaxation, the re- 
laxation path towards the minimum of the excited state 
PES (F^E) follows the minimum energy path, i.e. the 
dashed line drawn across points F and E in Fig. Ic. So, in 
spite of being driven by orientational degrees of freedom 
of the solvent, the dynamics of the system occurs along 
an effective coordinate that describes a concerted motion 
of solvation degrees of freedom and internal vibrations. 
Once again, the coupling between the two kinds of motion 
originates from their common coupling to the electronic 
system, and is a direct consequence of the molecular po- 
larizability. 

Time-resolved fluorescence spectra can be calculated 
following the same procedure described in the previous 
section for steady-state spectra, provided the temporal 
evolution of the Qo probability distribution is known. To 
such an aim we rely on a Fokker-Planck description of 
the dynamics of the joint probability distribution in the 
space of the solvation coordinate and momentum22ii£, 
W{Qo, Pq), as described by the following equation of mo- 
tion: 



dW A 9W ■ dW 



dW 



dP? r 



) 
. (7) 

where 7 is the phenomenological friction coefficient for 
the motion along the relevant relaxation path. Start- 
ing from point F, where the distribution is rigidly trans- 
lated from the equilibrium ground state distribution, the 
Fokker-Planck equation can be easily integrated along 
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the F-E path in Fig. Ic. The calculation can be simpli- 
fied by recognizing the overdamped nature of the solva- 
tion coordinate (7 ^ loq). In the overdamped regime, 
the momentum rapidly attains the equilibrium, so that 
it needs not to be considered as an independent dynamic 
variable^. The only relevant dynamics involves the Qq- 
distribution, w{Qq), as described by Smoluchowski equa- 
tion: 



dw 
'dt 



d^V dV dw 



kT 



dQl. 



(8) 



where V{Qq) is the (classical) potential for the Qo mo- 
tion. We verified that for 7 > 2ujq the Smoluchowski 
equation leads to the same results as the Fokker-Planck 
equation. Apart from a simplification in the integration 
procedure, working in the overdamped limit has the ad- 
vantage of reducing the number of free parameters re- 
quired to model solvation dynamics. In fact, as it is well 
knowK^" and as we have explicitly verified, in that limit 
all dynamic properties depend on 7 and ujq only through 
the ratio r — t represents an effective solvation 

time, and can be fixed as the longitudinal relaxation time 
(tl) of the pure solvenli. 

Fig. 3 shows time-resolved spectra calculated for DCM 
dissolved in two different solvents (CH3CN and CHCI3). 
All molecular parameters, and the solvent relaxation 
energies as well, are kept fixed at the values previ- 
ously obtained from the fit of steady-state spectra^!, 
(cf Fig. 2, left panels). According to literature data, 
we fix T = 0.3 and 2.8 ps for CH3CN and CHCI3, 
respectivelyii^. The spectra in Fig. 3 compare nicely 
with experimental spectraiii^, and accurately reproduce 
the progressive red-shift of the band, as well as its nar- 
rowing. In order to better understand the origin of the 
evolution of band-shapes during solvation. Fig. 4 shows 
time-resolved spectra calculated for the same parameters 
as Fig. 3A, but without accounting for inhomogeneous 
broadening, i.e. by collapsing the w((5o) distribution into 
a (5-function. The comparison between inhomogencously 
broadened spectra in Fig. 3 A and homogeneous spectra 
in Fig. 4 demonstrates that the evolution of the band- 
shape accompanying the solvent relaxation is essentially 
due to the evolution of the vibronic structure (i.e. of 
HR factors) , and reflects quite naturally the evolution of 
the molecular geometry along the relaxation path. The 
w{Qq) distribution, obtained from the solution of the 
Smoluchowski equation, also evolves with time, but it 
affects calculated spectra in a more subtle way. In fact 
w{Qo) is responsible for inhomogeneous spectral broad- 
ening, whose temporal evolution is hardly recognized in 
the fairly broad emission spectra of polar chromophores 
in polar solvents. 

Similar features, including the progressive red-shift of 
the emission band and its narrowing, are recognized in 
time-resolved emission spectra of many other dyes, in- 
cluding neutral dyes like coumarin 153^^, and zwitteri- 
onic dyes like LDS-750i2,. In fact, in spite of the oppo- 
site solvatochromism, neutral and zwitterionic dyes are 



predicted to have similar behavior with respect to time- 
resolved emission, as demonstrated by the time-resolved 
spectra calculated in Fig. 5 for the same zwitterionic dye 
whose steady-state absorption and emission spectra are 
reported in Fig. 2, right panels. Both the red-shift and 
the narrowing of time-resolved emission bands are due to 
the increase of the mixing between the two basis states 
during the relaxation of the system along the effective 
solvation coordinate, as discussed in the previous section 
with reference to steady-state spectra. 

Other time-resolved experiments can be described 
within the proposed approach. Transient hole-burning 
experiments have been designed to follow the relaxation 
along the ground state PES. Specifically, we make ref- 
erence to the experiment described in Ref. In this 
experiment, a solution of coumarin 102 (C102) dissolved 
in CII3CN is impulsively excited to transfer population 
from the minimum of the ground state PES (point G in 
Fig. lb) to the vertical excited state (point G' in Fig. Ic). 
After a delay time long enough as to allow for the com- 
plete relaxation of the system in the excited PES (from 
G' to E in Fig. Ic), a second dump pulse is applied to 
vertically transfer the equilibrated population in the ex- 
cited state PES into the ground state PES (from E to E' 
in Fig. lb). We are not interested in modeling the ini- 
tial relaxation of vibrational coordinates occurring dur- 
ing the first few tens of fs^, but, much as it was done for 
time-resolved fiuorescence, we want to model the slower 
spectral evolution due to the solvent relaxation. We then 
assume an instantaneous (within the relevant solvation 
timescale) relaxation along the vibrational coordinate, 
towards point R in the ground state PES (Fig. lb) lying 
on the local minimum path with respect to Qi. The cal- 
culation of transient absorption spectra then goes via the 
solution of the Smoluchowski equation for w{Qo) from 
time zero (where the equilibrium w{Qo) relevant to the 
excited state PES is rigidly translated to point R in the 
ground state PES) towards the ground state equilibrium 
in G. 

Figure 6 shows transient differential absorption spectra 
calculated for C102 in CII3CN. Molecular parameters for 
C102, as well as the solvent relaxation energy relevant to 
C102 in CH3CN, are obtained from the fit of steady-state 
absorption and emission spectra^. The additional pa- 
rameter entering the Smoluchowski equation, t, is again 
fixed to 0.3 ps. In order to compare directly with ex- 
perimental results, the upper and lower panel in Fig. 6 
report the spectra calculated in two different temporal 
windows. The comparison with experimental spectra in 
Fig. 5 of Ref. is very good: in particular we repro- 
duce not only the blue-shift of the band, but also the 
appearance of a TIP, and its temporal evolution. The 
origin of TIPs is easily recognized in the variation of the 
band-shape of the transient absorption that accompanies 
its progressive blue-shift. Whereas the rigid translation 
to the blue of the transient absorption band is quite eas- 
ily understood based on the simplest linear model for 
solvation, the temporal evolution of the band-shape is 
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a direct consequence of the molecular polarizability. In 
fact, as long as the system relaxes along the R-G path in 
Fig. lb, the molecule readjusts its polarity (from p ^ 0.23 
to 0.16, in the specific case of C102 in CH3CN), and HR 
factors relevant to the absorption process change with 
time. TIPs are puzzling features of time-resolved spectra 
that cannot be explained easily within the standard pic- 
ture for solvation, but are a natural consequence of the 
polarizability of the solute molecule. 

TIPs appear quite commonly in pump-probe 
experimentsi2i2242ii2^. In these experiments, after 
the application of a pump pulse that drives the system 
to the excited state PES, the variation of the optical 
density of the solution is monitored by a probe-beam as 
a function of the delay time after the pulse. The exper- 
iment is easier than the hole-burning experiment, but 
the interpretation of resulting spectra is more complex. 
In fact at least three processes can contribute to the 
observed signal at each frequency. Both the bleaching 
of the absorption from the ground state (due to the 
reduced ground state population) and stimulated emis- 
sion decrease the optical density, leading to a negative 
signal, whereas the absorption from the excited state 
increase the optical density, originating a positive signal. 
A reasonable separation of the different contributions is 
only possible if the three processes occur in well sepa- 
rated spectral regions. For polar chromophores in polar 
solvents, the large Stokes-shifts guarantee in general a 
reasonable discrimination between signals originating 
from the bleaching of the ground state absorption and 
from stimulated emission. If the positive signal from 
excited state absorption does not interfere too much, 
the pump-probe experiment gives information on the 
evolution of the emission band during the relaxation of 
the system along the F-E path in the excited state PES 
(Fig. Ic), so that, in favorable cases, the pump- probe 
experiment gives similar information as time-resolved 
fluorescence. According to the previous discussion, we 
then predict a progressive red-shift and narrowing of the 
emission band, in quite good agreement with available 
experimental data for several chromophores, with either 
a neutraliSiSiiS^ or a zwitterionic ground stateSS. We 
underline once more that the red-shift of the stimulated 
emission band is easily rationalized in the standard 
solvation model. On the opposite, the narrowing of the 
band, being related to the variation of the molecular 
polarity with the solvent relaxation, can only be under- 
stood in models where the molecular polarizability is 
properly accounted for. The concomitant red-shift and 
narrowing of the emission band are responsible for the 
appearance of TIPs in observed spectra. 



IV. DISCUSSION 

Push-pull chromophores are interesting molecules 
in several respects: they find application as solva- 
tion probesSi^, as laser-dyes, and as chromophores for 



non-linear optics—. They are also good two-photon 
absorbersiSr and work as molecular rectifiersii. From a 
theoretical point of view, the possibility to describe the 
electronic structure of these molecules in terms of an ef- 
fective two-state model makes these systems fascinating 
model systems to understand electron-transfer in solution 
and specifically to investigate the physics of solvation in 
electron-transfer processes^i^^*^Si^. In a series of recent 
pg^pgj,g9j25j27j28^ wc have extended the standard two-state 
model for push-pull chromophores^Si^2iffi to account for 
solvation and molecular vibration. The model allows for 
a non-perturbative treatment of the relevant couplings 
and fully accounts for the solute (hyper)polarizability, 
yielding to a successful picture for steady-state spectra 
of push-pull chromophores in solutionSSiSSi^. Here we 
extend the model to time-resolved spectra and demon- 
strate the important role of the solute hyperpolarizabil- 
ity in several characteristic and non-trivial spectroscopic 
features. 

Consider a simple system with no electron-vibration 
coupling. Due to the molecular (hyper)polarizability, the 
reaction field affects the solute properties and, in turn, is 
affected by the presence of the solute. As a consequence 
of this feedback interaction, the ground and excited state 
PES for the Qo motion are deformed^^. The linear po- 
larizability {ac/E) is responsible for a renormalization of 
the relevant harmonic frequencies^^: 



Vta/E = t^oy 1 - '^'^QaGjEl (9) 

while higher order polarizabilities originate anharmonic 
correction a?^i?? that can be fairly large for chromophores 
with large NLO responses. The role of the molecular 
first polarizability, a, on solvation dynamics is easily 
understood and has already been discussed by several 
authors^iiSSiSi^'*. If the hyperpolarizability is neglected, 
the PES are harmonic and the solvation relaxation func- 
tion, C(t), has a simple exponential form. The relax- 
ation time, TQjE — iI^'g/e^ however renormalized 
with respect to the relaxation time of the pure solvent, 
r = 7/^0 1 ^"^^ depends on specific properties of the so- 
lute. A qualitatively different behavior is expected if hy- 
perpolarizabilities are accounted for, leading to anhar- 
monic PES. In this case, a non-exponential C{t) is ex- 
pected, whose detailed behavior depends the details of 
the relevant PES. The role of the solute hyperpolariz- 
abilities on solvation dynamics has not been recognized 
so far and the observed deviations of C{t) from the sim- 
ple exponential behavior have always been ascribed to the 
failure of the continuum model for the solvent or to the 
underdamped nature of the solvation coordinateiii^^. 
Instead a non-exponential behavior can show up within 
the simplest model for solvation, as a consequence of the 
solute hyperpolarizability. 

The molecular (hyper)polarizability has even more im- 
portant effects in the presence of coupled internal vibra- 
tions. In fact the common coupling of vibrational and 
solvation degrees of freedom to the electronic system in- 
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duces a mixing between the two kinds of coordinates, 
even if no direct coupling between the modes is included 
in the hamiltonian. The physical origin of the mixing 
lies in the molecular (hyper)polarizability: the polariz- 
able molecule in fact readjusts its polarity to any vari- 
ation of solvation or of vibrational coordinates, but this 
variation of polarity affects in turn all coordinates, lead- 
ing to an interdependence of all motions. As discussed 
in the previous section, even in the simplifying hypoth- 
esis of a complete separation of the timescales of vibra- 
tional and solvation motion, the relaxation of the sol- 
vent coordinate always implies, for a polarizable solute, 
a concomitant relaxation along the vibrational coordi- 
nates. Therefore time-resolved experiments performed 
on polarizable chromophores measure a complex motion 
where the molecular geometry readjusts following the sol- 
vent relaxation. Deviations of the observed relaxation 
time from the relaxation time of the pure solvent, as well 
as deviations from the simple exponential behavior for 
C{t), are obviously predicted in this model. Just as an 
example. Fig. 7A compares the spectral correlation func- 
tion calculated for hole-burning experiment on CI 02 de- 
scribed in the previous section. Deviations of the exact 
result from the simple exponential behavior, as well as an 
overall slowing down of the solvation dynamics with re- 
spect to the pure solvent, can be appreciated from the fig- 
ure, even if, due to the large experimental uncertainties in 
C {t) , they are possibly too small to be safely identified in 
experimental data. Larger effects are predicted for more 
polarizable solute molecules, as shown in Figure 7B for 
phenol blue, an interesting solvato chromic probe whose 
large hyperpolarizability is known since date^. Quite 
irrespective of their magnitude, these effects are impor- 
tant since they demonstrate that the presence of a po- 
larizable solute affects solvation dynamics: the relaxation 
time obtained from the spectral relaxation function is not 
a property of the pure solvent, but also depends on the 
solute molecule. Good solvation probes should be largely 
solvatochromic and hardly polarizable as well. 

The two-state model unavoidably assigns the excited 
state a linear polarizability equal in magnitude and oppo- 
site in sign to the ground state polarizability, then point- 
ing to a solvation time in the excited state smaller than 
the longitudinal relaxation time of the pure solvent. As 
discussed in Ref. |5^ the negative polarizability of the ex- 
cited state is probably an artifact of the two-state model, 
that neglects the mixing of the low-lying excited state 
with higher excited states whose (positive) contribution 
to aE usually over-compensate the negative contribution 
from the mixing with the ground state. Estimates of 
excited state polarizabilities are available for just a few 
chromophores and all confirm positive values^2*S. 
However, the possibility of observing a negative a for the 
excited state of highly polarizable chromophores cannot 
be excluded a priori. In any case, it is possible to modify 
the proposed approach to account for the role of (a few) 
higher excited states. However, this would introduce new 
adjustable parameters in the model, whose definition is 



arbitrary in the absence of accurate experimental infor- 
mation on the polarizability of the excited state, or on 
the relaxation dynamics in the excited state. 

The analysis of C{t) is hindered by large experi- 
mental uncertainties, and the effects of the solute (hy- 
per)polarizability are more easily recognized in the tem- 
poral evolution of band-shapes. As recently suggested 
by Matyushov^, if one accounts for the solute linear po- 
larizability, the different curvatures in the ground and 
excited state PES lead to different equilibrium w{Qq) 
distributions in the two states. Then a temporal evo- 
lution of w(Qo) is easily predicted in time-resolved ex- 
periments, resulting in an evolution of the inhomoge- 
neous broadening profiles in observed spectra. The so- 
lute hyperpolarizability makes the phenomenon more 
complex: anharmonic PES, as originated by non- linear 
polarizabilities22i2^, in fact lead to non-Gaussian distri- 
bution profiles. However, it is important to recognize 
that, as discussed in the previous section, band-shapes 
are actually governed by vibronic profiles, inhomoge- 
neous broadening being only responsible for their smear- 
ing out. The observed evolution of band-shapes is in- 
deed originated by the coupling between solvation and 
vibrational degrees of freedom, as induced by the molec- 
ular (hyper)polarizability. The role of vibronic profiles in 
defining observed band-shapes was recognized in Ref. Isol 
but the temporal evolution of HR factors was introduced 
ad hoc. We demonstrate instead that it is a natural con- 
sequence of the solute (hyper)polarizability. 

The evolution of band-shapes that accompanies the 
red- (blue-) shifts of time-resolved fluorescence (absorp- 
tion) bands can be responsible for a crossing of ob- 
served spectra. The resulting TIPs are not true isos- 
bestic points, as their location smoothly varies with time. 
In any case, the appearance of TIPs has puzzled exper- 
imentalists and theorists since datei^iSSiSiS^. Recently, 
a paper appeared aimed to demonstrate the failure of 
the two-state picture for push-pull chromophores, based 
on the appearance of TIPs in pump-probe spectra^^. In 
this paper, TIPs are ascribed to a barrier-activated re- 
action occurring in the excited state towards a transient 
emissive product. Whereas further theoretical and exper- 
imental work is needed to settle this point, we underline 
that the molecules investigated in Ref. 22^ are among the 
best chromophores for NLO, and are therefore largely 
(hyper)polarizable. The observed large deviations of re- 
laxation times from the longitudinal relaxation times of 
the pure solvents, point quite nicely to deviations from 
the standard linear perturbative picture of the solute- 
solvent interaction and perfectly fit within our two-state 
model for polar-polarizable chromophores, and suggest 
an alternative explanation of TIPs. 



V. CONCLUSIONS 

Time-resolved spectra of polar-polarizable chro- 
mophores in solution are described based on a contin- 
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uum solvation model. A two-state picture is adopted 
for the solute, and linear coupling to internal vibra- 
tions is also accounted for. The model, originally de- 
veloped to account for NLO responses^i^ and steady- 
state spectreiSSiSi of push-pull chromophores in solution, 
is extended here to describe time-resolved spectra. The 
comparison with experimental data confirms that the 
proposed model quite successfully grasps the essential 
physics of polar-polarizable chromophores in solution. 
The model presents some obvious drawbacks. The two- 
state picture must be extended to account for higher- 
energy states if detailed information is required on ex- 
cited state properties. Moreover, the Holstein model 
for electron-vibration coupling only accounts for linear 
terms: quadratic coupling is expected to be important 
in cases when a very different geometry (and hence very 
different bond orders) can be attributed to the two ba- 
sis states. Again the extension of the model in this di- 
rection is easy, but requires detailed spectroscopic in- 
formation to fix all the model parameters. Finally, the 
continuum model for solvation has serious limitations, 
particularly in cases were large site-specific interactions 
(like, e.g., H-bonds) are expected. The hypothesis of 
complete separation of vibrational and solvation time- 
scales can fail for molecules with particularly slow inter- 
nal degrees of freedom and/or for solvents with a very 



fast dynamics. In spite of these limitations, and quite 
irrespective of its applicability to some specific systems, 
the proposed model opens the way to obtain exact, non- 
perturbative solutions of the solvation problem also in 
the presence of electron- vibration coupling. The spectral 
consequences of the molecular (hyper)polarizability, that 
is unavoidably neglected in the standard linear perturba- 
tive approach to solvation, are important and non-trivial, 
and show up quite consistently in steady-state and time- 
resolved spectra of push-pull chromophores. In partic- 
ular here we have demonstrated that deviations from 
the exponential decay of the spectral relaxation func- 
tion, that are usually taken as proof of the failure of 
the continuum model for solvation, appear quite natu- 
rally within the simplest picture, provided the molecular 
(hyper)polarizability is accounted for. Similarly, the tem- 
poral evolution of band-shapes, as well as the appearance 
of TIPs, usually ascribed to the failure of the two-state 
model, are another consequence of non-linearity. More 
complex models can of course be developed, and will lead 
to a better or more detailed description of spectral prop- 
erties of push-pull chromophores in solution. Based on 
an Occam-razor approach, we believe however that the 
available body of experimental data does not require any 
more complex physics than described here. 
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FIG. 1: (a) Potential energy surfaces calculated for the same parameters used below to fit spectra of DCM in CH3CN: 
zo = 1-14 eV, V2t = 0.88 eV, ei = 0.45 eV, wi = 0.16 eV, eo = 0.75 eV. Qo and Q\ describe the solvation and vibrational 
coordinate, respectively. The labels ahs and em mark the vertical lines along which the absorption and emission processes occur. 
For the sake of clarity, the energy gap between the two surfaces has boon enlarged, (b) Isopotential lines for the ground state 
and (c) the excited state PES shown in panel (a). Lines are relevant to equispaced energy values, the spacing corresponding 
to wi. The dotted lines drawn across the isopotential lines mark the equilibrium Qi value as a function of Qo- Letters on the 
graphs mark points discussed in the text. All axes are in arbitrary units. 
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FIG. 2: Left panels: (A) Absorption and (B) fluorescence spectra calculated for DCM (same parameters of Fig. 1) in hexane 
(continuous lines, eo ~ 0), CHCI3 (dashed lines, eo = 0.32 cV), CH3CN (dotted lines, eo = 0.85 eV). Right panels: (C) 
Absorption and (D) fluorescence spectra calculated for a zwitterionic chromophore {zo = 0.1 eV, V2t = 1 eV, ei = 0.5 eV, 
uji = 0.18 eV) for £0 = (continuous lines), £0 = 0.36 eV (dashed lines) and £0 = 0.85 eV (dotted lines). The intrinsic line- width 
for the Gaussian associated to each vibronic feature is cr = 0.06 eV. 




FIG. 3: Time-resolved fluorescence spectra calculated for DCM in (A) CH3CN and (B) CHCI3. The parameters are the same 
as for Fig. 2, left panels; t = 0.3 and 2.8 ps in panels A and B, respectively. 
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FIG. 4; The same as Fig. 3A, but without accounting for inhomogeneous broadening of the spectra. Again times relevant to 
the different spectra are marked on the graph. 
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FIG. 6: Transient diflterential absorption spectra calculated for C102 in CH3CN: zo = 1.28 eV, \/2t = 1.2 eV, ei = 0.33 eV, 
uji = 0.16 eV, eo = 0.6 eV. Spectra are relevant to the times reported on the graphs. Calculated spectra are shown in two 
separate windows for a better comparison with data in Ref. i22 . Arrows indicate the direction of increasing time. 
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FIG. 7: Solvation correlation function for the hole-burning experiment, calculated for (A) C102 in CH3CN (the same parameters 
as in Fig. 6) and (b) phenol blue in CH3CN {zo = 0.7 eV, V2t ^ 1 eV, ei = 0.42 eV, uji = 0.2 eV, eo = 0.7 eV, from Ref. 29). 
Continuous lines refer to the exact results; dashed lines correspond to the non-polarizable-solute approximation; dot-dashed lines 
have been calculated by taking into account the linear polarizability only; dotted lines correspond to the best single-exponential 
fit of the exact curve. 



